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Abstract. We deal with the viscous profiles for a class of mixed hyperbolic- 
parabolic systems. We focus, in particular, on the case of the compressible 
Navier Stokes equation in one space variable written in Eulerian coordinates. 
We describe the link between these profiles and a singular ordinary differential 
equation in the form 

Here V £ R d and the function F takes values into M. d and is smooth. The real 
valued function f is as well regular: the equation is singular in the sense that 
C(V) can attain the value 0. 



We focus on mixed hyperbolic-parabolic systems in one space variable in the 
form 

(2) E(u)u t + A(u, u x )u x = B(u)u xx . 

Here the function u takes values in $L N and depends on the two scalar variables t 
and x. We focus on the case the matrix B is singular, namely its rank is strictly 
smaller than N . In particular, a conservative system 

ut + f{u) x = (^B(u)u x ^j 
can be written in the form ([2]). Indeed, one can set 

A(u, u x ) = Df(u)-B(u) x , 

where D f denotes the jacobian matrix of /. In the following we will consider 
explicitly the case of the compressible Navier Stokes equation in one space variable: 

Pt + {pv) x = 
(pv) t + (^pv 2 + pj = (vv x ) 



(3) 



pe + P - 



1 ., 
-pv + pe + p 



= (kO x + vvVx^j 



Here the unknowns pit, x), v(t, x) and 9(t, x) are the density of the fluid, the 
velocity of the particles in the fluid and the absolute temperature respectively. The 
function p — p(p, 6) > is the pressure and satisfies p p > 0, while e is the internal 
energy. In the following we will focus on the case of a polytropic gas, so that e 
satisfies 

fAS R() 

4 e = 7, 
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where R is the universal gas constant and 7 > 1 is a constant specific of the gas. 
Finally, by v(p) > and k(p) > we denote the viscosity and the heat conduction 
coefficients respectively. 

In the following, we assume that system ([2]) satisfies a set of hypotheses intro- 
duced by Kawashima and Shizuta in [8]. These conditions were defined relying on 
examples with a physical meaning, in particular, up to a change in the dependent 
variables, they are satisfied by the equations of the hydrodynamics. For complete- 
ness, we recall these conditions here. First, we assume that the rank of the matrix 
B(u) is constant and we denote it by r. Also, B(u) admits the block decomposition 

( 5 ) = ( b (u) ) • 

The block b(u) belongs to M rxr and there exists a constant q, > such that for 
every |* E W 

(6) (&(«)£ > Q,|£?. 

In the previous expression, (•, •} denotes the standard scalar product in K. r . 

Also, we assume that for every u the matrix A(u, 0) is symmetric. We denote 

by 

A n (u) AUu) \ F( s _ ( E 11 {u) E^(u) \ 
A 21 (u) A 22 (u, u x ) ) £J[U) \ E 2l {u) E 22 (u) J 

the block decomposition of A and E corresponding to ([5]) . Note that only the block 
A 22 can depend on u x . Finally, we assume that for every u, the matrix E(u) is 
symmetric and positive definite. 

In the following, we focus on two classes of solutions of ((2|) : travelling waves and 
steady solutions. A travelling wave is a one variable function U (y) satisfying 

(8) [A(U, U') - aE(U)] U' = B(U)U". 

In the previous expression, a is a real number and is the speed of the wave. From 
(|5J| one obtains a solution of ([2]) by setting u(t, x) = U(x — at). Steady solutions 
are solutions of ((2]) that do not depend on time: namely, they are one variable 
functions U(x) satisfying 

(9) A(U, U')U' = B{U)U". 

We also require that U is bounded on x £ [0, +oo[ and admits a limit as x — > +00. 
Steady solutions in this class are sometimes called boundary layers. In the appli- 
cations, it is often interesting to focus on the case the speed a in ((SJ) is close to an 
eigenvalue of E~ 1 A. Since in general is not an eigenvalue of E~ 1 A : it is useful to 
distinguish between and (HJ). 

Travelling waves and steady solutions are powerful tools to study the parabolic 
approximation of hyperbolic problems: since the literature concerning this topic is 
extremely big, we just refer to the books by Dafermos [5] and by Serre [11] and 
to the rich bibliography contained therein. Concerning the analysis of the viscous 
profiles, we refer to Benzoni-Gavage, Rousset, Serre and Zumbrun [2], to Zumbrun 
[13] and to the references therein. Loosely speaking, the problem of the parabolic 
approximation of an hyperbolic system is the following: consider the family of 
systems 

(10) E(u £ )uf + A(u e , eu%)u% = eB(u s )u e xx , 



(7) A(u, u x ) = 
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which reduces to (f2|) via the change of variable u(t, x) — u e (et, ex). Letting e — > + , 
equation (fTOj) formally reduces to 

(11) E(u)u t + A(u, 0)u x = 0. 

As we pointed out before, the compressible Navier Stokes equation in one space 
variable can be written in a form like (|10[): in this case, system (|11|) . which formally 
is obtained by setting e = 0, is the Euler equation. The proof of the convergence of 
u £ to a solution of (fTTj) is an open problem in the most general case, but convergence 
results are available in more specific situations (see e.g. Bianchini and Bressan [3], 
Ancona and Bianchini [JJ and Gisclon [7]). Travelling waves have been exploited 
to study the parabolic approximation (fTD) and its relations with the limit (fTTj) . 
especially to describe such phenomena as the formation of singularities. The study 
of steady solutions has provided useful information to the analysis of (fTD)) especially 
in the case of initial boundary value problems. In the following we will be mainly 
concerned with viscous profiles having sufficiently small total variation. 

In a previous work [4] , the authors studied the approximation (|10p for an initial 
boundary value problem with special data. Concerning the structure of the matrices 
E, A and B, we introduced a new condition of block linear degeneracy, which is 
the following. Let the block An (it) and -Eii(m) be as in ([7]). The block linear 
degeneracy precribes that for every real number a the dimension of the kernel of 
[An(u) — aEn(u)] is constant with respect to u. In other words, the condition 
of block linear degeneracy says that such a dimension can in general vary when a 
varies, but it cannot vary when u varies. 

If the condition of block linear degeneracy is violated, then one might face 
"pathological" behaviors, in the following sense. In [4] Section 2] we discuss an 
example of a system that satisfies all the Kawashima Shizuta hypotheses, but does 
not satisfy the block linear degeneracy. We exhibit a steady solution ([9]) of this 
system which is not continuously differentiable. Note that in this case it is not 
completely clear, a priori, what we mean by solution of ([5]), because we are dealing 
with non regular functions. The details are given in [4]. 

On the other side, imposing the block linear degeneracy is restrictive in view 
of some applications. In particular, as pointed out by Frederic Rousset in [9], 
this condition is satisfied by the compressible Navier Stokes equation written in 
Lagrangian coordinates, but is violated by the same equation written in Eulerian 
coordinates. 

The problem is the following. Consider the Navier Stokes equation written in 
Lagrangian coordinates: 

(n - v x = 
( e+ y) t + H = (^ + ^),- 

Here the unknowns are r(t, x), v(t, x) and 6(t, x): r is the specific volume, and 
it is the inverse of the density, r = 1/p. As in (O, v and 9 denote the velocity 
and the absolute temperature of the fluid. As pointed out for example in Rousset 
[TO] , the compressible Navier Stokes equation written in Lagrangian coordinates 
can be written in a form like (f2|), with E, A and B satisfying all the Kawashima 
Shizuta [8] conditions. In particular, the rank of the viscosity matrix is constantly 
equal to two. Let An be the same block as in (fjj. In the case of the compressible 
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Navier Stokes written in Lagrangian coordinates An is a real valued function and 
it is actually identically equal to 0. This implies, in particular, that the condition 
of block linear degeneracy is satisfied by the compressible Navier Stokes written 
in Lagrangian coordinates. Indeed, let En(u) as in ([7]): in this case, En(u) is a 
strictly positive, real valued function. It follows that, for every real number a, the 
dimension of the kernel of — <jEn(u) does not depend on u. 

By direct computations one can then verify that the viscous profiles of the Navier 
Stokes equation written in Lagrangian coordinates satisfy an ordinary differential 
equation which does not have the singularity exhibited by |T]). 

Let us consider now ([3]), the Navier Stokes equation written in Eulerian coordi- 
nates. We want to write it in the form 

(13) E{u)u t + A(u, u x )u x = B(u)u XXl 

requiring that the matrix A(u, 0) is symmetric. In the following, we assume that 
p is bounded away from 0, say p > c p > for a suitable constant c p . This implies 
that the system does not reach the vacuum. 

We proceed as in Kawashima and Shizuta [8j and by multiplying ((3]) on the left 
by a suitable nonsingular matrix we eventually obtain that ([3|) can be written in 
the form (fT3"|) for 
















V 











k/9 



J 0\pe-w x /6 pvee/6 - k! 'p x /9 



b = I ( n uia ) an=P P /{0p)- 



(14) Elf, r. f)) = -\ I) i> U ) /J(/>. r. 

and 

I I P P v/p P P 

(15) A(p, v, 9, p x , v x , 9 X ) = - p p pv - v 1 p x p g 

\ pg - vv x j9 pveg/9 - k'p x /9 

In the previous expression, we denote by p p and pg the partial derivative of p with 
respect to p and 9 respectively, while by exploiting ((4]) we get that eg = R/{^ — 1). 
Also, to simplify the notations we write 

(16) A 21 = l(^) A^ = \ ( PV ~ V ' PX 
and 

k/9 

The condition of block linear degeneracy is violated here. To see this, let us focus 
on the case a — 0: the dimension of the kernel of An — a\\V is if v ^ 0, but it is 
1 when v = (we recall that p p > 0). 

To see what in principle can go wrong, we focus on steady solutions, the situation 
for travelling waves being analogous. We set 

(17) w = p x z= (v x , 9 x ^j . 

Then (J9j) becames 

a n v A^(u) \/w\/0 \ / w x 
A 2 i(u) A 22 (u, u x ) )\ z ) I b{u) JU 
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i.e. 



auvw + A 21 z ' = 
A 2 iw + A 22 z = bz x 



Assume d^O, then 



(18) 



can be written as 

A T 7 



auv 



= b 



-l 



A 



22 



A 2 iA^ 



auv 



Note that the matrix b is invertible and hence the previous expression is well defined. 

By combining (fT7|) and (p~8|) , one obtains that the steady solutions of the Navier 
Stokes written in Eulerian coordinates satisfy a singular ordinary differential equa- 
tion in the form 



(19) 

provided that U 



dU 
dx 



1 



F(U) 



F(U) 



T 

z I and 



V 



A^z/an 

v z 

A 22 v - A 2 iA 21 /au 



\ 
) 



We say that equation (|19| is singular since in general v, the velocity of the fluid, 
can attain the value 0. 

Let us now go back to the example in [4] we mentioned before: the example 
deals with a system with non continuously differentiable steady solutions. It turns 
out that, as a consequence of the fact that the system violates the block linear 
degeneracy, the steady solutions V satisfy a singular ODE in the form 

(20) d Z = 77^(n 



dt 



where C is a real-valued, smooth function that can attain the value 0. The loss of 
regularity experienced by V is actually due to the fact that PD|) is singular. 

Summing up, we have the following: the condition of block linear degeneracy is 
satisfied by the compressible Navier Stokes equation written in Lagrangian coordi- 
nates, but it is violated by the same equation written in Eulerian coordinates. As 
a consequence, the viscous profiles of the equation written in Eulerian coordinates 
satisfy a singular ODE and hence might in principle experience a loss of regularity. 
This suggests that we should look for more general conditions than the block linear 
degeneracy. Namely, we want to find a condition on the singular equation (|20|) 
which is sufficiently weak to be satisfied by the Navier Stokes equation written in 
Eulerian coordinates and by other systems violating the block linear degeneracy. 
On the other side, this condition should be sufficiently strong to rule out any loss 
of regularity. 

The definition of this condition and the corresponding analysis of the singular 
equation (|20|) is done in [5] . The conditions defined there apply to a class of systems 
that do not satisfy the block linear degeneracy and, in particular, to the Navier 
Stokes in Eulerian coordinates. 
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Remark 1. For a different approach to the analysis of the viscous profiles of the 
Navier Stokes equation in Eulcrian coordinates we refer for example to Wagner |12j 
and to the references therein. 
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